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Abstract 

In this paper we prove the Local Asymptotic Mixed Normality 
(LAMN) property for the statistical model given by the observation of 
local means of a diffusion process X . Our data are given by X s+i d/i(s) 
for i — 0, . . . , n— 1 and the unknown parameter appears in the diffusion 
coefficient of the process X only. Although the data are nor Markovian 
neither Gaussian we can write down, with help of Malliavin calculus, 
an explicit expression for the log-likelihood of the model, and then 
study the asymptotic expansion. We actually find that the asymptotic 
information of this model is the same one as for a usual discrete sam- 
pling of X. 

[Frangais] Dans ce papier nous demontrons la propriete LAMN pour 
le modele statistique constitue par I'observation des moyennes locales 
d'une diffusion X. Nos donnees sont definies comme X s+i dfi(s) 
avec i = 0, . . . , n — 1 et le parametre inconnu apparait seulement dans 
le coefficient de diffusion du processus X. Bien que cette observation 
ne soit ni Gaussienne ni Markovienne nous pouvons, par le calcul de 
Malliavin, obtenir une expression pour la log-vraisemblance du modele. 
Nous sommes alors capables de calculer I'information asymptotique et 
montrons qu'elle est la meme que pour I'observation ponctuelle de la 
diffusion. 
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1 Statement of the problem and main results 



1.1 Introduction 

Model. Let us consider the family of strong solutions to the fol- 
lowing scalar equation 



where -B is a one-dimensional Brownian motion. Wc suppose that 9 lies in 
some compact interval G of R and that is a real constant, which does not 
depend on 9 and thus is known to the statistician. 

Observations. We consider ^ some probability measure on [0, 1] and 
assume that our observation of the process is given by the local means of X 
associated with this measure, with sampling of size 1/n: 



(observations) Xj = Xj^n-= / Xs+j dju(s), for j = 0, . . . , n — 1. 

Jo " 

In the sequel this case is referred to as the integrated diffusion case. This 
is an indirect observation of the process X and the observation is no more 
the realization of a Markov chain. Thus, this framework is deeply related 
to the inference of hidden processes. We assume that /x does not depend 
on 6 and is known by the statistician. When /x is equal to the Lcbcsgue 
measure, the observation is the discrete sampling of It = /q Xgds. This 
is presumably the simplest case of the observation of only one component 
of a bidimensional diffusion process {Xt, It)o<t<i-, which is known in the 
literature as the standard integrated diffusion case. Clearly, the usual case of 
pointwise observation of X is obtained if fi is some Dirac measure. However 
we will exclude that the measure has mass only on the end points of the 
interval and hence make the assumption: 



This paper is concerned with the Local Asymptotic Mixed Normality prop- 
erty of this statistical model. 

Motivation. Taking as the observation the integrated process is actu- 
ally quite natural. For instance, it arises when the realization of the process 
has been observed after passage through an electronic filter. Also, in random 
mechanics (see Kree and Soize [19]), X models the velocity of the system 
and in general, we observe its position, i.e. the integral of X. The modeling 
of ice-core data can be made through an integrated diffusion process (see 
Ditlevsen, Ditlevsen and Andersen [2]). Integrated processes also play an 
important role in finance, when modeling the stochastic volatility (see for 
instance Barndorff-Nielsen and Shephard [1] and references therein). 



dX^ = a{X^,9)dBt + b{X^)dt, 



(1) 
(2) 





m((o,i)) >0. 



(3) 
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Literature background. Despite of these numerous motivations, few 
statistical studies deal with this situation. Gloter [7] [8] provides an esti- 
mator in the multiphcative case a{x, 9) = 9a{x) and proves its consistency 
and asymptotic normahty. The case of a low frequency observation (local 
means over interval of length 1) is studied by Ditlevsen and S0rcnsen [3], 
using prediction-based estimating functions. On the other hand, for a direct 
observation of the diffusion X, there are many contributions in the liter- 
ature: see Genon-Catalot and Jacod [5], Prakasa-Rao [24] and references 
therein. None of these works deal with the problem of optimal estimation 
in the integrated diffusion model. 

Here, we directly address the problem of the LAMN property, whose fun- 
damental consequence is to provide information on the minimal dispersion 
for an estimator of the parameter (see Ibragimov and Has'minskii [14], 
Jeganathan [16] [17], Prakasa-Rao [24], Le Cam and Lo Yang [21]). Such 
properties, for the observation of a discrete sampling of the diffusion, have 
been established in the one-dimensional setting by Dohnal [4] , and then ex- 
tended by Gobet [10] [11] to the multidimensional setting, both in the high 
frequency and ergodic framework. For this, Malliavin calculus techniques 
were used and paved the way to possibly handle more general situations than 
Markovian observations. This is exactly this way we follow in this work, to 
tackle the case of integrated diffusion. 

Outlook. We guess that this model captures the main difficulty of most 
hidden models: the lack of Markov property for the observation. Hence 
the method developed below (augmented observation, Malliavin calculus 
representation, Gaussian approximation) may be useful to treat more general 
situations. Among the natural situations coming from applications, one 
can think of the measurement of a stochastic phenomenon blurred by some 
noise, or stochastic volatility models widely used in finance [6] . This can be 
formalized as follows: the system is governed by the d + d'-dimensional 
stochastic differential equation 

Jo Jo 

where only a discrete sampling of the first d components is observable. This 
is left to further research. 

1.2 Main results 

Before going into the details of our results, we present a very simple example 
which gives some insight on the type of results that one can expect. 

Example 1 (Multiplicative Brownian case). Assume that the model is 

Xf = OBt 

(corresponding to b = and a(-,6) = 9,^o = 0)- 
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1. Consider a first situation where one observes the diffusion at discrete 
times. Hence, the observation is {Xi/n)o<i<n> or equinalently [Zi = 
6{Bi/ji — B(^i_iyji) = OGi)i^i<n! where Gi are independent centered 
Gaussian variables, with a known variance. Thus, the estimation of 
9^ is achieved at rate -s/n, with a minimal variance equal to 20'^. 

2. Now consider a second situation where one observes only the inte- 
grated diffusion at discrete times. Hence, the observation is (Xj = 
6 Jq B (s+i) fi{ds) = OG'j)i<i<n, where {G[)i is a centered Gaussian vec- 

n 

tor, with a known covariance matrix. In addition, this matrix is invert- 
ible and thus, 6^ can be estimated with the same rate and asymptotic 
variance as before. 

This means that observing the process at discrete times or its integrated 
version lead to the same accuracy in the parameter estimation. The results 
of this paper state that this is true, even for the more general models (l)-(2), 
which is far from intuitive. 

Before stating our main results, wc define the working assumptions of 
this paper. The coefficients a : R x ^ R and h : M ^ M, are assumed 
to satisfy the fohowing set of conditions (as usual, derivatives w.r.t. are 
denoted with a dot: for instance, dga = a). 
Assumption (R) 

1) The function a : R x 6 ^ R is C^+t for some 7 G (0, 1) (it admits 
a derivative which is 7-H61der). The one dimensional functions x 1— > 
a{x, 6), X i—i- d{x, 9), X !—>■ b{x) are assumed to be C^(R). 

2) The functions a, a and b and all their derivatives with respect to x are 
bounded uniformly in 9. 

3) We have the non degeneracy condition, for some a: a{x, 9) > a > 
for all X, 9. 

Actually, the uniform controls in (R) can be weakened to local ones, using 
extra techniques of space localization (see Lemma 4.1 in [10]). Wc omit 
further details. An extension of our results to a multidimensional parameter 
9 and to time dependent coefficients is straightforward, in the same way as 
it is done in [10] and [11]. 

We denote by the law on C([0, 1]) of the process , and then simply 
denote X the canonical process on C([0, 1]). We let p"'^ denote the law 
on R" of the observation O" := (Aj)j=o,...,n-i, when the true value of the 
parameter is 9. And for two values of the parameter we introduce the 

likelihood ratio, 

dp"'^i 
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The main result is that this statistical model satisfies the so called LAMN 
property. For this denote the sequence ?i„ := n^^/^, and let Oq E Q and 
/i G M such that + ^n^i G ©j Vn. Then, by the following theorem, the 
model has the LAMN property for the likelihood at point Oq, with rate u„ 
and conditional information: 

Theorem 1. Assume (R), then we have the expansion, 
log Zl^eo+ur^h = - 1/2/124 + R^, 

where In > , Rn and there exists an extra random variable 

n — >o<j n — >OG 

N ~ A/'(0, 1) independent of the process X such that, Nn converges in law 
under to Ny%^. 

Moreover this convergence is stable: for any random variable F measur- 
able with respect to the X, we have the convergence in law (F, N^) ' 

n— »oo 

(F, N ^Tqq). In particular it implies the joint convergence under P^°; 

n— >oo 

Remark 1. Let us stress that the rate Un = n"^/^ and the information Iq^^ 
are the same one as for the pointwise observation (see Genon-Catalot and 
Jacod [5]). This corroborates the intuition from Example 1. 

We will not be able to prove directly this result, instead we shall consider 
first the easier problem where one can observe additionally the exact value 
of the diffusion at some instants. This device was proved to be useful in 
Gloter and Jacod [9] for the study of a Gaussian diffusion process observed 
with noise that leads to non Markovian observations too. 

Let = fen be an integer in {1, ... , n} and define L = L„ := [n/fej , then 
we consider the set of random variables: 

(^n,aug ^ O'^u |Xm,/ = U{Xi}. 

Since this set of variables contains more data than the initial set, we call 
it the augmented observation set. Clearly, we can split this set into blocks. 
Bo, . . . ,IBl, where for / = 0, . . . , L - 1 

= {Xki, ■ ■ ■ ,Xki+k-i,Xm^iyn} 

and = {XkL, ■ ■ ■ ,Xn-i,Xi]. Note that if kL = n we consider that the 
last block is empty, and (immediate) modifications should take care of this 
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in the sequel, however to have shorter notations we will not explicitly write 
these modifications. 

The advantage of this set of augmented observation is that using the 
Markov property of X, the law the block B; conditional to the previous 
blocks (B;/)/'</ only depends on the last variable, Xh, of the block B/_i. 

n 

Denote by p^'^^S'^ the law of O^-^iig on M"+^+i and introduce the like- 
lihood ratio for the augmented observation: 

j„n,aug,ei 

^n,aug ^ op .^n,aug) _ (5) 

Theorem 2. There exists a sequence A;„ — ^ oo, such that the augmented 
model satisfies the LA MN property: 

log = HNr - i/2h'ir + Rr, 



where , ^ and there exists an extra random vari- 

able N ~ A^(0, 1) independent of the process X such that, Nn^^ converges 
in law under to N -^/Tq^. Moreover this convergence is stable. 

^From Theorem 2 and from the consequences of the LAMN property, an 
asymptotically optimal estimator On in the augmented model should satisfy 
that y/n{9n — 6Q) is asymptotically distributed under Pg,, as —^=N. However 

any estimator in the initial model of observation O" can be seen as an 
estimator in the augmented model, hence the Theorem 2 is sufficient by 
itself to imply a lower bound for estimation in the initial model. 

Remark 2. The fact that fc^ — oo means that the data added in the ob- 
servation are sparse compared to the initial data. Actually, the Theorem 2 
holds for any sequence kn whose growth to oo is slow enough. 

If wc assume now that = G N remains fixed as n ^ oo, the number 
of data {Xki)i=o l added to the model is not negligible compared to the 

n 

number of initial data. Hence the statistical properties of the augmented 
model shall depend on k and thus differ from the statistical properties of the 
initial model given in Theorem 1. Actually we have the following LAMN 
property for the augmented model in that case. 

Theorem 3. // a sequence kn = k is fixed, then the augmented model sat- 
isfies the LAMN property with rate Un = n~^l'^ and conditional information 
equal to: 



{X„eo)ds. 
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As expected, the conditional information is greater by a factor {k + l)/k, 
due to the non neghgibihty of the added observations. Actually this factor 
should be read as 1 + ^, meaning that an addition of ^% of data increases 
the information in the same way. Local means and values at discrete points 
are not redundant (as expected from the multiplicative Brownian case, see 
Example 1) and moreover, they bring an equal information. Considering 
A; = 1 is interesting, since we observe then on each block [i/n, (i + l)/n] 
both the exact value Xj^ and a mean Xi. It appears that the asymptotic 
information is then twice the information given by the observation of only 
the exact values (or only the means). 

1.3 Outline of the paper 

In Section 2 we study the score function given by the observation of only one 
block of data (Bq for instance). We first focus on the existence of a density 
for a block of data; and in the case of a block of size 2, (n^/^ jj^ (X^/^ — 
Xo)djLt(s),n^/^(Xj_ —Xq)) we give original lower and upper bounds of Gaus- 
sian type for the density. It is useful for our proof of the LAMN property, 
but it is also interesting for itself. 

In Section 2.2 we present an exact expression for the score function of 
a block of data Bq (see Theorem 5). This result is the key point in the 
proof of the LAMN property, it extends a former result of Gobct [10] [11] 
which gave the score function for the observation of X\ . In Section 2.3 we 

study an explicit approximation for the score function when the sampling 
interval 1/n tends to zero and the length of the block k/n remains moderate 
so that one can consider the coefficients of the diffusion X almost constant 
on the interval [0,fc/n]. The key point is the Gaussian approximation for 
the diffusion given in Section 2.3.1. 

In Section 3 we deduce from the previous results a proof of Theorems 
2-3 and Section 4 deals on how to deduce Theorem 1 from Theorem 2. 

Finally the Appendix contains the proof of some results of Section 2.1 
together with some useful lemmas. 

Notations. In our proofs, we will keep the same notation for constants 
which may change from one line to another. In particular, the constants 
c,c{k),c{p),c{p,k) will stand for all finite, non-negative and non-decreasing 
deterministic functions of an index p (arising from L^-norm) and of the block 
size k. These constants are independent of n, 9 and depend on the process 
X^, only through the bounds on the coefficients a, b and their derivatives. 
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2 Score function for a block of data 



In this section we shall study the law of the blocks of data Mi; recalling 
the Markov property of the process X it is sufficient to focus on Bq = 
{Xq, . . . , Xk-i,Xk/n} assuming that the diffusion X now starts from some 
value Xq. In this section it is convenient to transform the short time asymp- 
totic k/n ^ into an almost stationarity property of the coefficients. To 
this end, we introduce the rescaled process X"' = n^{X\ — xq) (where X^ 

n 

solves (1) with Xq = xq). It solves the equation 

dX,^'' = aniXl"'', e)dWt + bn{X,^'')dt, = 0, (6) 

where is a standard Brownian motion (arising from the rescaling of B), 
and 

a„(x, 9) = a{xo + n"^/^x, 6*), hn{x) = n~^/%{xQ + rT^/'^x). (7) 

Since for the score we are only concerned with the law of X'^'^ , we can 
assume that W is independent of the rescaling coefficient n. 

2.1 The density of an integrated diffusion 

In this section, we will present preliminary results on the density of the law 
of the mean of a diffusion process. However the proofs are postponed to 
Section 5.1. To our knowledge, the lower and upper bounds for this density 
are new results. 



2.1.1 Existence of the density 

Our first result actually deals with the two dimensional variable given by 
solely one local mean and the exact value: 

^^n,e^ yn,e^ . ^ X^''dii{s),xf'^ (8) 

fnV2 /"'(x^, -x,)d^i{s),n'l\x\-x,)\ . 
\ Jo " / 

Notice that, by the Markov property, the preliminary study of this bi- 
dimensional variable will be a key step to obtain results on the observation 
vector O". 

Theorem 4. Assume (R), then the vector (C/"'^,y"'^) admits a density 
Pxq{-, on MP, and there exist two constants ci > C2 > 0, such that, 

The constants ci and C2 only depend on the bounds on the coefficients a, b 
and their derivatives. 
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The proof of this theorem is given in Section 5.1. The existence of the 
density is obtained by means of the MaUiavin calculus. On the other hand, 
the upper and lower bounds rely on the direct study of {U^'^ , V"''^) around 
its skeleton (see Hirsch and Song [12] [13] for related works; and Kohatsu- 
Higa [18] for different methods involving Malliavin calculus). 

The following is a direct corollary of Theorem 4: 

Corollary 1. The vector = (Xq, . . . , X^-i, admits a positive 
density. 

Proof. The bi-dimensional process (-'^z, );=o,...,ifc-i is a Markov chain 

n ^ 1 

with transition density {xi,xi+i,6) = np^^ (na (xi+i-xi), {xi+i-xi),d). 
Then it is clear that the vector Bq admits a positive density. □ 



2.1.2 Invertibility of the Malliavin covariance matrix of a block 

Actually the existence of a density for the law of the random variable Bq will 
not be sufficient, and we need a non degeneracy condition for this variable. 
Before this, let us precise briefly a few notations from the Malliavin 

calculus, used in the sequel (see Nualart [22] [23] for details). Wc let H 
be the Hilbcrt space L^([0, oo)) so that the Brownian motion {Wt)te[o,oo)j 
appearing in (6) , is canonically associated to this Hilbert space via the stan- 
dard isometry. In this setting, for any p > 1 and natural number q, 
recall that the set B*?'^ denotes the space of real valued Wiener function- 
als with q derivatives and whose derivatives belong to L*'(J1). If wc de- 
note by D the derivative operator then the space O'^''' is endowed with the 



norm, H-f' 



mm + Y.UE(\\D'F\\l,t,,.,,) ' . The space of vari- 



L2{[0,oo)i) 

able with q derivatives in any L^'(O) is denoted W'°° = np>iD^'^. These 
deflnitions can be extended to random variables with values in any Hilbert 
space V and the corresponding spaces are denoted W'P{V), W'°°{V) (see 
Section 1.5 in Nualart [22]). In particular the operator D is then well defined 
from D'''°° to B)^^^-'^(H). Finally, the adjoint operator of D is the Skohorod 
integral S, and the Malliavin covariance matrix of an element F G D^'°° 
is defined as the matrix jFi,-,Fa = [{D-Fi, D.Fj) 
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Now, we consider the variables, 



K''--= t K^'Ms) (10) 

[/r'^= [\Ki'i-K'')Ms) (11) 

Jo 



K\ ■■= - KiU)Ms) (12) 

^0 



- K;U)Ms). (13) 

^0 



Note that the joint law of these k+1 variables is, by rescaling, the same as the 
law of the vector composed with variables of the first block Bq: n2(XQ — 
xo,Xi — Xq, . . . ,xl._i — X^_2,Xfe — xl._i). These variables satisfy the 

n 

following non degeneracy property whose proof is postponed to Section 5.1.3. 

Proposition 1. Under (R), {Uq'\ . . . ,U^''^) Denote by K{e) the 

Malliavin covariance matrix of {Uq'^ , . . . , U^'^)- It is a.s. an invertible ma- 
trix and for all p > 1, we have 

E{\detiKi9))\-P)<cip,k). 

2.2 An exact expression using Malliavin calculus 

In this Section we intend to give an exact expression for the score function 
of the observation of Bq or equivalently for the vector {Uq'^ , . . . , U^'^) given 
by (10)-(13). 

Under the condition (R), we know that there exists a version of the 
solution of (6) such that P— almost surely the function 6 is continu- 

ously differentiable for all t and r"'^ := ^'^g is a solution of the stochastic 
equation (see Kunita [20]): 

^^n,e ^ ^(A-f e)TfdWt + ^(A-f '^ e)dWt + ^{X^^')T^''dt (14) 

To"'^ = 0. 

The main result of this section is an explicit representation for the derivative 
of the log-likelihood of one block. This extends a former result given by 
Gobet (see [10] [11]). 

Theorem 5. The random vector {Uq'^ , ■ ■ ■ , U^'^) admits a positive density 
on Mf''^^, denoted by Pxoiuo, ■ ■ ■ ,Uk,&)- For a.e. {uq, . . . , Uk), this density 
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is an absolutely continuous function with respect to the parameter 9 and we 
have the formula: 



PXQ 

Pxo 



{uo,...,Uk,9) = E 



vO<ij'<fe 



de 



(Uj' = ■Uj)j=0,...,A 



where K{9) ^ is the inverse of the Malliavin covariance matrix of {U^'^ , ■ ■ ■ ,Uf 

Proof. Denote f/"'^ the Wiener funetional, t/"'^ = {Uq'^ , . . . , Uj!'^) and let 
/ : — >• M be a smooth function with compact support. Then the 
function 9 ^ E \_f{U^'^)\ can be differentiated pointwise and: 



Tnfi\ 



89^ 



f{U 



E 



n,0 



j=0 



89 



By Proposition 1 the Malliavin covariance matrix of C/"'^ is invertible and a 
standard computation on Wiener functionals (see formula (2.4) p. 81 in Nu- 

alart [22]) shows that: |J(f/"'^) = T!}'=o{D{f{U^'')),DU;/)^K{9)^^.,. 
It follows that ^E [/(C/"'^)] is equal to 



E 



BIT 

Y^Y.{Dif{U-'')),DU;'')^Ki9)j],. ^ 

3=0 j'=0 

where is the if- valued random variable 



n,0 



89 



E 



D{f{U^'')),L' 



H. 



k k f)TT'"'fi 



K{9)-}DUp' 



j=o j'=0 

Introducing 6 the adjoint operator of D, we get 

8 



89 



E 



E 



f{U"''')6{r 



(15) 



(16) 



Let g be any smooth function with compact support on M. Using the inte- 
gration by part formula and the equation (16) we have: 



J d9g{9)E{f{U'''')) 



d9g{9)^E [f{U-'') 



■ J d9g{9)E [f{U^'')6{LC 



d9gi9)E [/(C/-'^)S[5(L'^) | (U^'' , [/,"''')] 
Introducing the density of the random vector U"''^ the equation above writes, 

J g{0)d9 J f{uo, . . . ,Uk)pxoiuo, ■ ■ ■ ,Uk,9)duo . . .duk 
= -J 9{d)d9 J f{uo, . . . ,Uk)E[S{L^) \ {Ul"'^ = ui)i]pxo{uo, ■ ■ ■ ,Uk,9)duo . . .duk- 
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Now using Fubini's theorem it can be seen that duo . . . dw^-almost 
everywhere the function 6 — >■ Pxo{uo, . . . ,Uk,0) is absolutely continuous with 



Pxo{uo,---,Uk,d) = E[S{L^) I (?7"'^ = ui)i]pxoiuo,...,uk,0). 
Hence the theorem is proved. 



□ 



Remark that the proof of Theorem 5 does not rely on the specific expres- 
sion (10)-(13) and thus an analogous representation for the score function 
seems achievable in many situations. 

2.3 A Gaussian approximation for the log-likelihood 

In this section we intend to give a tractable approximation for the score 
function of (U^'^ , . . . , U^'^)- 

2.3.1 Approximation for the diffusion 

We introduce = a{xo, 9)Wt and ff = a{xo, 6)Wt which stand -by (6) and 

(14)- for the first order approximations of -Y"'^ and r"'^ = ■ Then, we 
consider the quantities obtained by replacing in (10)-(13) the process X by 
this Gaussian approximation 

:= a{xo, e) f Wsd^i{s) = a{xo, 9) C ^i{[s, l])dW^„ (17) 

JQ Jo 

Ul := a(xo, e) !\wj+s - VFj-i+.)dM(s), for j = 1, . . . , - 1 (18) 
Jo 

= a{xo, 9) J^' ^ fi{[0, s-{j- l)])dWs + a{xo, 9) j'^ ^l{[s - j, l])dWs, 

:= o(a;o, 9) [\wk - Wk-i+s)dfi{s), (19) 
Jo 

= a{xo,9) [ fi{[0,s-{k-l)])dWs, 

Jk-l 

where we have repeatedly used the Fubini theorem for stochastic integrals 
(see [25] p. 176). In the next lemma we control the difference between the 
Uj ' and their approximation in terms of Sobolev norm. 

Lemma 1. For all k,p > 1, there exist constants c{k,p), c{p) such that for 
all j G {0, 



dul 

89 



U 



~d9' 



2,P 



< c{k,p)n 



-1/2 



3,p 



< c{p), 



< c{k,p)n-^^^, 



2,P 



~d9 



< cip), 



(20) 
(21) 



3,p 



12 



VO<i,j'<fc, 



E ( Uf'uY - U^UP 



< c{k)n 



-1 



(22) 



Proof. The inequalities on the right hand side of (20)-(21) are immediate 
by the definition of Uj. 

Comparing expressions of (10)-(13) with (17)-(19), the two remaining 
bounds in (20)-(21) will be a consequence of the Minkowski inequality - for 
the Sobolev norm - and of the control on the diffusions: 



sup 

t<k 



+ sup 

2,P t<k 



n - t: 



2,P 



We only prove the control on X"''^ since the proof for r"'^ is analogous. 
Recalling (6)-(7), we can write 



^>h9 _ p^9 



1 [' ^ uXl 



\an{X':-^ e) - a{xo, e)]dWs + f 6„(A';'^)ds 

^0 



,e)X^^'^AudWs + -^ I bixo + '-^)ds. (23) 

But we know [22] that under (R) the variables X^*-'^ belong to D^'^ with a 
control (independent of ^,n): sup„^ „2<s<fc ^(I^mi m2'^"' I^) — "^(P'^)- This 

< n~^^'^c{p, k) after a few computa- 



is sufficient to deduce 
tions. 



2,P 



To obtain (22) note that by (20) it is sufficient to show E {^{Uj'^ - U^)Up^ < 

c{k)n~^. This property will follow again from an analogous relation on the 
diffusion, 

sup E ({Xl'''^ - X^)X^,) <c{k)n-^. 

t,t'<k ^ ' 

Indeed, from (23), the above expectation is equal to 



n 



i-tAt' fl 

-V2 / E 
Jo Jo 



E[a'^{xo,e)X^'^] = J^a'^{xo,e)E[bn{X::''^)]du < cn'^/^ E [b{xo)Wt>] = 



I'^ixQ + n-^/'^uX^'^, e)X^'^] dua{xQ, eo)ds+ 

f E \h{xQ + n-^l'^X^'^)wA 0(^0, 0)ds. 



Using 

and the boundedness of a'^^ and 6', we get the required estimate. 



□ 
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2.3.2 Approximation for the log-likelihood 

Let us denote the deterministic tridiagonal matrix K of size {k + 1) x (k + l), 

vi c O' 
C Vl + V2 ' ' ■ 

■•• ■•• ■■• 



K 








■ • Vl +V2 C 
C V2 



where the entries of the matrix are: 



Vl 



= f\{[s,l]fds, V2= f\{[0,s])''ds, c= f\{[0,s])fi{[s,l])di 
Jo Jo Jo 



It can be easily checked that a'^{xo,9)K is the covariance matrix of the 
Gaussian vector (C/q , . . . , U^) and that it is invcrtiblc using (3). Now the 
idea is to introduce the score function that would be produced from the 
observation of this Gaussian vector. Hence we let: 



Cxoi^o,- ■■,Uk,0) = -{xo,e) < 



a{xo,6) ^ ^ UjK^ j,Uf - {k + l) 

0<j,j'<k 



}■ (24) 



In this section, we will show that this quantity is an approximation for the 

true score function ^: 
p 

Theorem 6. Let us consider the difference, 

— {uo,.. .,Uk,9) -Cxoiuo,-- ■,Uk,0) := r^oiuo, ■ ■ .,Uk,9). (25) 



Then we have the following bounds: 



E 



r,,{U^'\...,Uf,e) <c{k)n-\ 



Vp > 1, E 



r,M\...,Uf,e) <c{k,p)n-"\ 



(26) 
(27) 



Proof. Keeping in mind the definition of (see (15)), we introduce its 
approximation based on the Gaussian quantities defined above: 

3=0 j'=0 

The first step is to obtain the following control on the difference ri := 

Vp > 1, ||ri||^i,p(^) < cik,p)n-'/\ (28) 
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Actually, it is a easy consequence of Lemma 1, Proposition 1 and the in- 
vertibility of K, noting that the Malliavin covariance matrix of coincides 
with the covariance matrix a?{xQ, d)K of the Gaussian vector . We omit 
further details. ^ 
The second step is to obtain a simple expression for 5{L^). To see this, wc 
first use the relation for F G B^'°°,u G 0^^°°{H), d{Fu) = F5{u)-{D.F,u) jj 
(see [22]): 

^(^') = E E -^a{x^,e)-^KT}5{D{U'^,)) 



, d9 

3=0 j'=0 



j=Oj'=0 \ I H 



On the one hand, S{D{Uj,)) = Uj, {SoD is the identity operator on the first 

chaos space). On the other hand, one has -q^ = a{t"o,e) ^j 
deduce 

d{xo,6) ^ 



'^'''•'^j:j:a{.«,e)-%}{D.UlD.Ul,)^ 



a{xn, 9) 

j=0 j'=o 

^ E E Ufaixo, erKr^uS, - + 1). 



j=0 j'=0 

Now set 



and take the conditional expectation in the relation S{L^) = S{L^)+5{ri): by 
Theorem 5, we get (25) with r^^{uo, . . . ,Uk,e) = E (^5{ri) \ iUj'^)j = {uj)j^ + 

E{r,\{Uf)j = {uj)j). 

The final step in the proof is to show that r^n, satisfies conditions (26)- 
(27). For the first condition, since the Skorohod integral has zero mean, we 

have E[rxo i^o '^i • • • ' ^fc ^)] ~ E{r2) and we conclude using (22). 

We now prove (27) . The conditional expectation being a contraction on 
it is sufficient to prove 

E{\5{ri)\^)p < c{p, k)n-^/\ E{\r2n^ < c{p, k)n'^/^. 
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The first estimate follows from (28) and the continuity of the operator 5 
from I])^'P{H) to I/. The second one is an immediate consequence of Lemma 
1. □ 

Remark 3. Let us note that the constants c{k),c(k,p) in Theorem 6 should 
increase as the block length k goes to infinity since the Gaussian approxima- 
tion ceases to he valid in that case. However in the sequel we shall not need 
a precise evaluation of this dependence on k since we will have the possibility 
to conveniently choose the growth rate of k = kn- 

In the following sections we will need this corollary of Theorem 6. 

Corollary 2. We have for all p > 1, 

E ^-^{u^\...,u]:^e)']<c{k,p). 



Proof. By Theorem 6 it is sufficient to show that E 



c,,,{u^'',...,u;:'",e) 



c{k,p). But from the expression of Cxq, this estimate is clear. 



< 



3 Asymptotic study for the augmented model 

In this paragraph wc establish Theorem 2. Let us recall some notations: we 
now deal with the diffusion given by (l)-(2); kn is some integer in {1, ... , n}, 
Ln = [n/kn\ and our observation consists of the Ln + 1 blocks Bq, . . . ,Bl„ 
described in Section 1. The length of the block is kn^i + l, where kn,i = kn 
if / < Ln — 1 and kn^L^, = n — Lnkn- For sake of simplicity in the sequel we 
sometimes omit the dependence with respect to n and I of the block size, 
and let kn,,i = k with a slight abuse of notation in particular for the last 
block of data. 

To be able to use the results of the Section 2, we introduce on each block 
the random variables corresponding to the definitions (10)-(13) for the first 
block. Hence for Z G {0, . . . , we define the kn^i + 1 following variables: 

Uo,i = n^Xkl-Xki), 

n 

Ui,i =n^Xki+i- Xki), 

Uk-i,i = {Xki+k-i — Xki+k-2), 
Uk,i = na (X fc(i+i) - Xki+k-i)- 

n 

Clearly the observation of the {Uj^i) for Z G {0, . . . , Ln},j G {0, . . . , kn,,i} is 
equivalent to the observation of the L„ + 1 blocks. Using the Markov prop- 
erty for the process X it appears that the law of the vector (J7jy)j=o,...,fc„ , 
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conditionally to all the variables C/^y with I' < l,j G {0, . . . , is the 

same as conditionally to X,t//n only; moreover this law - conditionally to 
^ki/n = xq - coincides with that of the vector {Uq'^ , . . . , Uj^'^) studied in 
Section 2. Thus it admits the density pxj^iiuo, . . . ,Uk,0) studied in Sec- 

tions 2.2-2.3. Hence the log-likelihood of the augmented model admits the 
additive structure: 

i21 PXki{Uo,l,---,Uk,l,0O + Unh) 

'— ^ n 

feo+u„h pXkiiUo^i, ■ ■ ■ ,Uk,i,s) 
^7^0 Px^{Uo,i,...,Uk,i,s) 

Owing to Theorem 6, we deduce the decomposition 

pOo+Unh 



1=0 •'^o 

Jen " 



In the above decomposition, we will show in Sections 3.1-3.2 that the explicit 
term involving Cxq governs the asymptotic behavior of the log-likelihood 
ratio; the other term does not contribute in the limit. 

3.1 Proof of Theorem 2: the expUcit term 

Let us introduce a slight modification of Cxki > which has the advantage of 

"n" 

being a smoother function w.r.t. 9: 



ii,n{e) = -{XMM{ 

a n 



0<j,j'<k 



}, (29) 



and we set N^^ = n„ Ef=o ^U^o) and iT^ = -ul Ef="o %(^o). 
Proposition 2. If kn ^ oo slowly enough, 

1=0 



^ / Cx,, {Uo,i, Uk,i, s)ds = hNr - + Rn, (30) 



where 1^"^ ^ ° > Xe„ , Rn ^ ° > and there exists an extra, random variable 

n— >oo Ji^oc' 

~ AA(0, 1) independent of the process X such that, Nn"'^ converges stably 
in law under to N ^/XZ,. 
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Proof. Comparing (24) with the definition of ^i^ni^) above and using a Taylor 
expansion for ^i^n{Q) around we get the equation (30) with a remainder 
term i?„ = + satisfying: 



^n^=E / [-{Xu,s)--{Xu,eo) 
'—0 ^" 



a n 



^ a(Xu , s)2 



1=0 




(k + l) 
(31) 
(32) 



(2) 

(fori? n we have used that 9 i-^ d(x, 9) is 7-H61der continuous). To complete 
the proof, we repeatedly use the following classical convergence result about 
triangular arrays of random variables. 

Lemma 2 (Genon-Catalot and Jacod [5], Lemma 9). Let {xi')o<i<L„, U 
he random variables, with being J^l'j^^-measurable. The two following 

conditions imply Yli=o xf ~* U : 



E E [x?\Jf] and ^ ^ [(xD'l^] ^ 0. 

1=0 1=0 

• We first focus on Nn^^. Let us introduce the sigma field J^p = a{XQ; Bg, s 
^) for / = 0, ...,L„ and = (^iXo;Bs,s < 1). Then the variable 

^l,n{9o) is .F/Yi -measurable and the asymptotic behavior of Nn"^ will fol- 
low from Lemma 2. To make clearer this point we introduce the following 
approximation based on conditionally Gaussian variables: 



< 



Cl,n{9) = -{XH,eo){ 
a n 



0<j,j'<k 



}■ (33) 



Here, Uj^i is the Gaussian approximation under of Uj^i corresponding on 
the block Mi to the variables (17)-(19) on the block IBq: 



Uq^i := a{Xki,9)n2 / (B ki+s — Bm )du(s) 



Uj^i ■.= a{Xki,9)n^ f {B 
Jo 



ki+j+3 — B ki+j-i+s )diJ,(s) for j = 1, . . . , A; — 1, 



Uk,i ■= a{Xki,9)n2 / (Bkg+i) — B ki+k-i+s )du(s). 

n Jo n n 

Observe that this vector (t(jy)j=o,...,fe has, under and conditionally to 
Xki = xo, the same law as the vector (C/^),=o ... k defined in Section 2.3. 
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Thus its conditional law is Gaussian with covariance matrix a{Xki,9)'^K. 

Hence, the variable ^;,n(^o) is J^j^i-Taeasurahle and under P^", it is condi- 
tionally (to Xki) distributed as a recentered x^(A; + 1) variable. Thus we 

n 

deduce the following four properties: 



1) UnJZtoEeo 6,n(^0)|^ 



0; 



2) Using = ~ n/kn oc and oo, one has 



4E^^0 te,n(0o))'|^ =ulY,^{K + l)(-] (Xm,^0) 



1=0 



^^''''^^\l\l)\xs,9o)ds + oMl) (34) 

2 



3) uU2iZo^9o 6,n(6'o) 1^/ 
to oo slowly enough; 



>2 / - {Xs,eo)ds=l0,; 



< cn '^Lnkfi^ < cn ^k^ — > 0, if kn goes 



0. 



^Prom these four properties, it follows (see Jacod [15]) that Un ^f=Q 6,n(^o) 
converges stably under P^" to a mixed Gaussian variable as in the statement 
of the proposition. To obtain the limit for A^n"^, it is sufficient to prove that 



(35) 



Due to Lemma 2 a sufficient condition consists in the two following points: 
1) UnEtoEeo \^l,n{^o) - ^lA^o) \ Jf] ^ in probability; 



2) ulEtoEoo iCwM-^UWl^^r 



in probability. 



But these two points can be shown using (20) and (22) of Lemma 1 (for ky, 
slowly increasing). 

• We now study In^^. A direct differentiation of ^;,n(^) (recall (29)) gives 



o<jj'<fe 



a n 



Then, with a few computations similar to the study of N^^, we obtain (for 
appropriate kn)'- 
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1) ulY^Eoo 4n(^o)|^^ =ulY,-'^ikn,i+l)^{Xu,9o)+OM'-^) 



I, I 1 -iXs,do) ) ds + Opeo(l) 



(36) 



2) KEtoEeo [ti,nm?\^i 



< cn ^kt 



0. 



Combined with Lemma 2, these two convergences imply that of /^"^ to 
under P^o. 



The remainder term Rn- Firstly a direct use of (20) gives E{ 



c{kn)n 



-7/2 



o(2) 



< 



if kn slowly goes to oo. Secondly the convergence to zero of 



i?n ' = Z^^o ^nl more delicate and Lemma 2 is helpful for this. To this 
end we evaluate the conditional expectation of R^] using (22) and the fact 
the {Uj^i)j have the conditional covariance matrix a{Xki,0)^K: 



9o+unh ^ ^ a{Xki,6o)'^ 

[-{Xkl,s) {Xkl,9o)]{ " 



6o 



-l}ikn+l)ds 



+ 0(n ^■Unc(A;„)). 



The function a being C^'^'^ in 6*, one gets: J^b^o ^doi-^ni I -^i 



C(fcn) -1/2 



for appropriate kn- With similar considerations we evaluate 



the second conditional moment and obtain X^^^q ^Oo 

c(kn)LnV^'^'^ "^°°> 0. 



?(l)^2 



□ 



3.2 Proof of Theorem 2: the negligible terms 

It remains to prove that, as announced, there is convergence to zero of 



J2i=oVl with r]i = Jg° rXkiiUo,l,---,Uk,i,s)ds. We aim at applying 

n 

Lemma 2 by computing the first two conditional moments of rji under F^°. 
The main difficulty here comes from the fact that we do not have an explicit 
expression for rx^^{{uj)j,9). Indeed by Theorem 6 we know bounds for the 
moments E^^^ (kxo((f^i)i) ^)D where by E^^ q we denote the expectation 
with respect to the law of X'^'^ solution of (6) . This is a priori insufficient to 
compute the conditional moments of rji under which involve quantities 
such as Eq ^ {Vxo{{Uj)j.,s)\^^) for s / ^o- Thus in Lemmas 7-8 in the 
Appendix we study the transformation of such moments under change of 
measure. 
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Firstly, we evaluate the conditional expectation of rii, 

f'9o+u„h 

But ElJr,{{U,)„s)]\ < \ElArAmj,s)]\ + \E^^Jr,i{U,),,s)] - El[,[^^^^^ 
can be bounded using (26) and Lemma 8 in the Appendix by c{k)n^^ + 
\s — 6o\ E'g^^[\rx{iUj)j, s)\"']'^ for some a > 1. Then by (27) we deduce 
\E[rii I :F"]| < c{kn)[unn~^ + «^n~^/^]. Finally a block length kn slowly 

increasing guarantees X^^^q Eoq [vi I -^f ] ^ 0- 

Secondly and similarly, owing to Theorem 6 and Lemma 7 in the Ap- 
pendix, we get E [r]f \ J^p] < c{kn)u^n~^ — 0. Therefore, by Lemma 2, we 

have proved Yli=o Vl ^ 0- This ends the proof of Theorem 2. 

3.3 Proof of Theorem 3 

The proof is essentially the same as that of Theorem 2, the difference in 
the asymptotic information comes from the difference in the limit of the 
quantities (34) and (36) when k is fixed. 

4 LAMN property for the initial model 

In this Section we are back to the model where the observation is only 
O"' = (Xj)j=o,...,n-i and we will prove Theorem 1 by relying on the LAMN 
property for the augmented model. 

A first intermediate result is that one can approximate the log-likelihood 
of the augmented model by a function of the observation O"". 

Proposition 3. There exist random variables r„ measurable with respect 
to O" such that: 

Proof. We have seen in Section 3 that ln(Zgji'J'^^_^„^^(0"''^'^s)) = hN^^ - 
1/2/i^/n"^ + Opeo(l) where the quantities A^n"^ and /n"^ were defined in 
Section 3.L 

Thus the proof of the proposition consists in introducing a proper mod- 
ification of these quantities which only depends on the observations. We let 
foi I = 0, . . . , kn 



a-\Xki-i,e) Yl Ui,,kj}Ui^j, - {k 
i<j,j'<k~i 
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with the convention X_i = is the known initial value of the diffusion and 
the matrix a^(xo, 6)K is the covariance matrix of the conditionally Gaussian 
vector (t/f, ,C/|_,): 



K = 



Vl + V2 
C 















C V1+V2 



Clearly, only depends on the observation O" since we have sup- 

pressed all occurrences of the variables Uq^i and Uk,i and replaced Xm 

n 

by Xi^i_i in the expression of (,in{0) (compare with (29)). Then we let 

• Study of Nr^^ — N°i^^. The first step is to consider the conditionally re- 
centered chi square approximation of ^i^{0) that we define as: 



i<j,j'<k-i 



Ui,jK7},Uij> - {k 



(37) 



The first step is to prove the validity of the approximation: 



(38) 



This is done similarly to the proof of N^"^ — Un '^i'=QCl,n{^o) ^ in propo- 
sition 2, by considering the first two conditional moments, but here the first 
moment is more delicate to handle: the conditional moment [C°^'^(^o) ~ 



^tni^o) I is of the form {k-l){g{Xki] 



-g{Xki-i)}h{XM-i)+0{c{kn)/n) 



Xkl 

n 

obs / 



X 



kl-l 



LP 



for g and h two functions. If we abruptly use the relation 

c(p)n-V2 then we only deduce that ^x„Ef=o^eofe°n(^o) "C^^o) I re- 
mains bounded in probability. To show that it actually converges to zero, 
we have to apply again Lemma 2 to the new triangular array of variables, 
■"n Y4=Q{k-'^){g{Xki)-g{Xki-i)}h{Xki-{). Then by rather long computa- 



< 



tions, using that - Xfe/_2| 

1 



LP 



< c(p)(A:/n)V2 and Eg,[Xui - X^i-i \ 



cn , we can prove, 

Un^{k - 1) 

1=0 



Ee,[{giX,_i) - g{Xki-i)}hiXki-i) \ < c{k)n 



-1/2 



0, 



, Y,{k - lfEe,[{g{Xu) - g{Xki-i)fh{Xki-if \ < c{k)n-^ ^ 0. 



1=0 
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Thus we deduce Ef=o ^eo[C?o) " C'(^o) I -^rl ^ 0- The second 
condition ul Ef=o EeMH^SOo) - ^(^o))'' I ^] < c(fc)n-i ^ is easily 
obtained and we deduce (38). 

Thus, in view of the equation (35), it remains to prove that n„ X^^oi'^?n''(^o)~ 
^/,n(^o)} is negUgible. But by Lemma 10 in the Appendix, comparing expres- 
sions (33) and (37) , it appears that conditionally to J^i' the random variable 
^irii^o) ~ Cl,n{%) is a recentered x^(2) variable and hence the following 
properties hold: 



1=0 

pn / (■ ^ ^ -.2 \ • 2 

^nE^^o({C(^o)-M^o)} \H) =^^X^4^(XM,eo)<f -0. 



These two properties imply by Lemma 2 the convergence to under of 

Ef=o {e/°n (^o) - e/,n(^o)}, and thus N^"- - iV°^^ ^ 0. 
• Study of IT^ - if"-. Exactly as we proved that tends to Xq, we can 
show that I'f^ Zqq. Thus the difference is negligible. 

Finally the proposition is obtained by setting r„ = hN!^^ — h? . □ 

Then Theorem 1 is a consequence of the following proposition combined 
with Proposition 3 and Theorem 2. 

Proposition 4. We have the convergence, 

Proof. The starting point is the relation between the likelihood of the ini- 



tial and of the augmented model: Z^^^g^^^^,^ = Eg„ [^^"'J'e^+u^ft I ^"J • By 

Proposition 3 we can write ^Q^^g^j^u^h ~ ^^"e'^" where En tends to zero in 
Pg(j probability. Using that r„ is measurable we deduce, 

Zl,eM - e^" = Ee, [e^-{e^- - 1) | O^] . 

We now use the inequality |e" — 1| < (1^1 A l)(e" -|- 1) to obtain that 
^Ifi^+ur^h - e^" I < an + with: 

an = Ee, [(|e„| A l)e^- \ O"] = Eg, A 1 | O"] e^-, 
(5n = Eg, [{\en\ A l)er»e- | 0«] = Eg, [{\en\ A 1)Z,-j7„%„, | 

It now remains to show the convergence to zero of a„ and 
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For an, let us notice that {e^")n is a tight sequence and that Eg^^ [|£„| A 1 | O"] 
converges in L^(Pgp) norm to zero since, 

Ee, [Eg, [\en\ A 1 I O"]] = Eg,[\en\ A 1] 0. 

For Pn, we have Eg,[Pn] = -E'6»o+«n/j[knl ^ But the sequence of 
probabihties P^" and P^o+«n/i restricted to the sigma fields O'*'^'^^ are con- 
tiguous (this is a consequence of the LAMN property for the augmented 
model, see e.g. Proposition 1 in Jeganathan [16]); hence the sequence (en)n 
which is measurable with respect to O"''^^^ and converges to zero in P^o- 
probability converges also in P^o+""'*-probability. This implies Egg[Pn] = 

E9o+Unh[\Sn\ A 1] ^ 0. □ 



5 Appendix 

5.1 Proof of results of Section 2.1 

Since the results of Section 2.1 concern only the study of a density for fixed 
values of 9, we omit the dependence upon 6 in our notations. We will prove 
the results in the following order. First in section 5.1.1, we show that the law 
of the Wiener functional (?7",y") = ^/J A'"d/i(,s), A'f ^ admits a density. 
Then we prove the lower and upper bounds given in Theorem 4 (section 
5.1.2) and eventually we deduce the Proposition 1 (section 5.1.3). 



5.1.1 Existence of the density p"^ 

We know [22] that under (R) the random variable -Y" is an element of ]D)^'°° 
and its first derivative is equal to 

where is the solution of 

dy^ = a'jxnyt<iWt + b'jxny?<it, 3^o" = l- (40) 
In the sequel we will repeatedly use the positivity of and the control 

Ei sup C^;,")^) + Ei sup (y^rn < cip). (41) 

tG[0,l] t€[0,l] 

^From this we can sec that the random variables U"' and are elements 
of ]D)^'°° and using (39) with the linearity of the operator D, we have 

AC/" = i{t<s}ynytr'an{xndf^{s) = a„{xn{yn-\t<i} [ y^d^is), 

Jo J[t,l] 

DtV" = anmynyn-\t<i}- 
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Using Theorem 2.1.2 p. 86 in [22], a sufficient condition for the existence of 

a density for (?7", F") is that its Malhavin covariancc matrix "yu'^y^ satisfies 
a non degeneracy condition given, for instance, by the foUowing lemma. 

Lemma 3. ^u^y^ is an a.s invertible matrix and for all p > 1, we have 

E{\detijur^yn)\-P) <cip). 

Proof. To have shorter notations, during the proof we will denote by c* any 
generic positive random variable which satisfies E{cZ^) < c{p). By direct 
computations we have, 



y:d,,{s)\dty^, 



[t.i] 

n\2 



Jo 

Now, define the probability density on [0, 1] 

< = <miy?r' (^J^imiy^r'ds 

and set /"(t) := Jj^ ^ 3^"d//(s). Thus we can write: 



(42) 
(43) 
(44) 

(45) 



deti^unyu) = {v\v^)l{y^)-' 



wllf^{tfdt 



rKr{s)ds 



Hence the above bracket can be interpreted as the variance of the function 
f"'{t) under the probability measure m"dt and hence: 



deti^u-y^) = {v^,v^)l{y^)-^ 



nt) 



m';f'\r)dr 



-I 2 



dt. 



But clearly under Assumption (R), (1^", F")^ > mft^[o,i]{yt)~^ infig[o,i](3^f )2 
and hence by (41) this yields, {V'^,V'''% > c^,, using our convention about 
generic positive random variables c*. Similarly, by (45), we have > 
and thus, 

2 

dt. 



dt, 



det(7t/",y") > ^ nt) - ^ m^r{r)dr 



Then, writing the integral above as 



.1/2 



mP/"(r)dr 



+ 



r(t + l/2)-^ m;?r(r)dr 
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and using the simple inequality x +y > {x — y) /2, we get: del{ju^yn) > 

/o^"^^ (/[t,t+ 1) 3^?d/x(s)^ dt. Using again inf^gp,!] > c*, we obtain: 

det(7[7n yn) > c* Jq^^ iJ,{[t,t + 1/2))'^ dt. But this integral is positive as soon 
as /i((0, 1)) > which is the case by assumption (3). Thus the lemma is 
proved. □ 

5.1.2 Bounds for the density 

For the proof of (9), we make a crucial use of the fact that the diffusion 
process -Y" is one dimensional by introducing the classical transformation: 

sn{x) := r a-\y)dy, := 
Jo 

By the assumptions on a, the function s„ is one to one on R and the 
derivatives of s„ and are bounded independently of n. By Ito's for- 
mula, solves the equation dW" = dWt + 6„(yV")dt where bn{w) : = 
^ o s~^{w) " ia'„ o s^i(w) and the initial value is Wq = Sn(A'^) = 0. We 
let P be the probability defined on ($7,^) by 



dP 
dP 



exp 6„(W:)dT^„ - i ^^(>V:)dn) 



The Girsanov theorem implies that the process W" is under P a standard 
Brownian motion. Note that the random variables (i/^jF") have the fol- 
lowing expressions with respect to this P-Brownian motion: 

/■\-i(W,")d/x(r), (46) 
Jo 

V- = s-\wn (47) 
Now let ho and be non negative real functions, then: 

Ep [/io(C/")/ii(F")] = Ep [/io(?7")/ii(y")L"] , (48) 

where = exp (^J^ 6„(W;?)d>V;? - ^ Jq bl{W^)dr^ . But using Ito's for- 
mula, = exp (Pn(>Vf ) - 5 Joibl + b'J{W?)dr^ where B„ is the primi- 
tive function of 6„ vanishing at zero. Since 6„ and b'^ are clearly bounded 

< 



by cn for some constant c only depending on a and b and 
cn"''-/^ |x| we have: c~-'-exp (— cn""-*-/^ |W"|) < < cexp 



Sn(x) 
Wi^l). By 

(47) and the boundedness of we deduce c"^ exp (— cn~^/^ |^"'|) < -^^ < 



^1/2 
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cexp we obtain: 



ho{U^)hi{V^)e 



n\ -cn~^\V\ 



<Ep[ho{U'')hi{y'')]<cEp 



/io(C/")/ii(^")e 



Hence we have transformed the problem of finding bounds for the density 
of the law of V'^) under P into an analogous problem under P. Conse- 
quently the bounds for p^^^ stated in (9) will follow from the next lemma. 

Lemma 4. Let ho, hi he two non negative functions. There exist some 
constants ci > C2 > 0, depending only on the coefficients a and b such that: 

c^^ J J ho{u)hi{v)e-''^^''''+''^Uudv < 

Ep[ho{U'')hi{V'')]<c^^ J J ho{u)hi{v)e-'''^'''+''^'>dudv. 

Proof. We first show the lower bound. Using that the random variable 
is measurable with respect to Wf (by (47)), we can write: 

Ep [ho{U^)hi{V'^)] = Ep [hiiV")Ep [ho{U^) I Wf]] 

giw)hi {s-\w)) Ep [/io(t/") I = w] dw, (49) 



where q is the density of the standard Gaussian law. Now let us admit 
temporarily the following relation on the conditional law of W^: 



Ep [^o(C/") I Wf] > c-^e-^W)' j ho{u)e-'''''du. 
Then Ep [ho{U'^)hi{V'')] is greater than: 

c-^ j ^o(^i)e~^"'dM X j Q{w)hi {s-\w)) e-'^'^^dw. 



(50) 



The change of variable v = s„^{w) in the second integral above, the inequal- 
ities \w\ < c\v\ and s^(u) > c give the new lower bound 

c-^ J /io(u)e-™'d« X y" {sniv)) hi {v) e'^^'^dv, 

with a new constant c. Since g is the Gaussian kernel and thanks to the in- 
equality |s„('u)| < c \v\, we deduce the required lower bound for Ep [ho{U"')hi{V"^)]. 

We obtain the upper bound quite similarly. Let us temporarily admit 
that for all e small enough there exists c{e) such that: 

Ep [/io(C/") I < c(e)-^e"W)' j /io(?x)e-"(^)"'du. (51) 
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Plugging this in equation (49), we deduce that Ep [ho{U^) \ Wf] is smaller 
than 

c(£)-^y /io(^i)e-'(")"'d« X J g{w)e''^'hi {s-\w))dw. 

Since q{w) = exp{—w^/2)/-\/2Tr, any choice of e smaller than 1/4 implies that 

the second integral in the equation above is bounded by c J e~^'^ hi (s^"^(k;)) du;. 
As for the lower bound, we conclude by the change of variable v = s~^{w). 

□ 

It remains to show (50)-(51). This is done in the following lemma. 
Lemma 5. For some constant c > and e > Q, we have 

Ep [hoiU'') I > c-^e-^'^^"^' J ho{u)e-^^' du. (52) 

For all £ g)0, e(, there exists c{e) > such that, 

Ep [hoiU'') I Wr] < c(£)-ie"(^r)' J /io(^i)e""("^"'du. (53) 

Proof. Let us recall that the process W^* := — tW" is a Brownian bridge 
on [0, 1], independent of the variable Wf . Thus, we can evaluate the condi- 
tional expectation Ep [/io(C/") | Wf = w] as the expectation (recall (46)), 



E 



ho I^J^s-' {W:+tw)dfi{t) 



(54) 



for W* some Brownian bridge. This Brownian bridge itself admits a decom- 
position 

W;=^m + Wr, (55) 
where ^ is a 7V(0, 1) variable, r] is the deterministic triangle shaped function: 



t if [0,1/2] 

(1-t) if [1/2,1]' 



and W** is the process on [0, 1] constructed as the concatenation of two 
independent Brownian bridges, one on [0, 1/2] and another on [1/2, 1]. Fur- 
thermore in this decomposition the r.v. ry and the process W** are indepen- 
dent. 

For any realization of W** we can introduce the real function, 

x^5VV**(^)= / {s-^ {xrjt + W;* +tw)}dn{t). 
Jo 

Using (55) and the independence of ^ and W**, the quantity (54) now writes, 

£;(W--)%)[/io(5W"(0)] (56) 
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where the inner expectation denotes the expectation with respect to the 
random variable ^ and the outer one with respect to the process W**. 

First we evaluate the inner expectation. Using that ^s a standard 
Gaussian variable we have 

E^^) [hoigw** im = (27r)-i/2 J hoig^,* (x))e-^ dx. (57) 

Note now that for any realization of W**, the function x i— > g'w**(a;) is 
differentiable and using that ^ < {s~^y < c we get 

- / r]tdii{t) < g'y^**{x) <c I ritdjiit). 

C JO JO 

By assumption (3) on the measure ii the integral ritdfi(t) is positive. Thus 
the function x <—>■ gyv** (x) is invertible on R, with a derivative bounded from 
above and from below by some constant independent of W** and n. This 
allows us to make a change of variable in (57) to obtain the bounds 

c / ho{u)e 2 du < [/io(5w"(0)] <c / ho{u)e 2 du. 

(58) 

Now the proofs of (52) and (53) are treated separately. 

• For the lower bound, we have seen that g^** is globally Lipschitz with 
a constant independent of W** and thus |5'^l,(^i)| < c|m| + |5yy',*(0)| < 
c\u\ + c\gy\>**(0)\. In addition, a simple computation from the definition of 

and then boundedness of (s~^)' show that 

|^**(0)| <c[|u;|+ sup iwri]. (59) 

te[o,i] 

Using this in (58) we find a new lower bound for the inner expectation: 

Ei^) [ho{gw**{0)] > c-^e-'^\-''^'''^W\^^r? j /to(u)e-'^"'dn. 

Taking the expectation with respect to W** proves that (56) is larger than 
c-iS(W") (^e-'^'^P^eidiW*)'^ X e-'="''//io(?^)e-™'dn. This gives (52). 

• For the upper bound, we write using that gy\!** is Lipschitz | l-ul — l^vv** (0)1 1 < 

\gw**{gy^**{y)) — fl'>V**(0)| < c|(7^^,,(n)|. Together with the relation (x — 
y)^ > x^^-ey^ (for x, y € M, £ € (0, 1)), we deduce that exp (-^(^^^(n))^) 
is upper bounded by 



where we have used (59). Combining this with (58) and taking the expec- 
tation with respect to W**, we get that the quantity (56) is smaller than: 

2 

c / ho{u)e~^S^e^'^^ Eyv** (^e'^^^P'eio.iiCWD^^ . The last expectation is finite 
as soon as £ is small enough, and thus (53) holds. □ 
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5.1.3 Proof of Proposition 1 

To have shorter notations we set ^ = /J A'J^jd)u(t), for j > 0. First we 
prove: 

Lemma 6. Let us define V'^ the Malliavin covariance matrix of the vector 

\ (X" , X^,^)] of size 2k. Then this matrix is a.s invertible and 

^ ■' J j=0,...,k-l 

E{det{T^)-P) < cip, k). 

Proof. In the case k = 1, the lemma reduces to Lemma 3. For k > 2, we 
proceed by induction by estabhshing simple relations between the columns 
of r*^ (this simplification follows from the flow property of the process A""). 

To see this, notice that firstly by (39) if t < k — 1 and s > k — 1, 
we have DtX^ = A'^'fc-i3^"(3'fc-i)~^ and secondly if t > fc - 1 and s < 
k — 1, DtX^ = 0. Using these two properties, a calculation shows that if 
(Cj)j=i^...^2fe denote the columns of T^, we have the relation, 



[C2k-1', C2k] 



J^' yii-i+siy]:-ir'Ms)^ c2k-2;y]:iy]:-ir'c2k-2 

/o 0\ 



+ 





V bk] J 



where 7^ is the matrix of size 2x2 given by 

\jl,{D,X-){DXk-i)dt jl,{D,X-fdt ) • 

This proves that det = dct Y^^^ dct 7^. But it can be seen that the matrix 
7/j has an expression similar to 7[/n,v'" (but with integration interval shifted 
from [0, 1] to [k — 1, k\) from which we can prove £'((det 7^)"*') < c{p). 
The lemma then follows from induction on A;. □ 

Now we can deduce the Proposition 1. Recalling (10)~(13) we can find 

an invertible matrix M of size 2k x 2k that maps ( A:"",, ) ) into 

V J^"^ J j=o,...,k-i 

a vector whose k + l first components are exactly {Uq, . . . , UJ^). Denoting 
the Malliavin covariance matrix of the image byMof (;f ^A^IYJ 

we have f^' = MT'^M*. Thus, Lemma 6 yields E({det{t'')yP) < c{p,k) 
since M is invertible. Observing that the Malliavin covariance matrix K{9) 
is the matrix extracted from the k+l first rows and columns of we deduce 
Proposition 1. 
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5.2 Some estimates on the change of measures 

For this section we denote by X the canonical process on C([0, oo)) and we 
consider the random variable on this space defined hy H = f(Uo, . . . , ^7^), 
where {Uq, . . . , Uk) is given by (10)-(13) with the canonical process X in 
place of X^'"'; we denote by Eg^^ the expectation with respect to the measure 
induced on the canonical space by the law of A'^'" solution of (6). 

Lemma 7. There exist r > 1 and a constant c{k) > 0, such that MH = 
f{Uo, ...,Uk)>0, ye, 9' e G, Vxo e M, we have 

E^',.,[H]<c{k)El,^ml. 

Proof. Recalling the notation of Section 2.1.1 we denote p"o(n, v, 0) the den- 
sity of the vector (8) and for j = 0, . . . , A; — 1 we let 



(60) 



Then using the Markov property of the process X under the laws Pq and 
Pq,, we have 



Xo 



k-l 
j=0 



< El,^[Hr^ El,^ 



j=0 



where r and r' are conjugate exponents. But we know by Theorem 4 that 
there exist two constants < C2 < ci (uniform w.r.t. 6,xo,n) such that 



Then one can bound the conditional expectation E^ ^^[{Z^^i q qiY' \ Xs,s < 
{k - 1)] by 



But if r is chosen large enough such that r' is sufficiently close to 1 the latter 
integral converges and is equal to some constant k. Proceeding by induction 
we get: 



E^ 



j=0 



^ Ee,xo 



j=0 



which gives the result. 



□ 
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Lemma 8. There exist c{k) > and a > 1 such that = f{Uo, . . . , 17^) 
(with E^^^JH]'' < +oo), ye,e' e e, Vxo e R, we have 



\E^,,,,[H] - El,^[H]\ < c{k) \e - e'\ [ElJHn^. 
Proof. Using the notations of Lemma 7, we write 



(61) 



k-l 



'k-i 

n ^ififi 



1 H 



j=i+i 



Thus for conjugate exponents a and the left hand side of (61) is bounded 

by 



k-l 



i=Q 

k-l 



k-l 



j=i+i 



i=0 



\Zi,e,0' - if E^^^^ 



k-l 



II {Zj,e,e'f \Xs,s<i + l 
j=i+i 



Using the Markov property of X it can be shown exactly as in Lemma 7 that 
the conditional expectation in the equation above is finite, as soon as /9 is 
small enough and bounded by k'^~^~^. Thus by Lemma 9 below, we deduce 

— J i k—i—1 



Eg,,,,[H]-Ee,:co[H] <ciP)\e-e'\ 



□ 



Lemma 9. There exists /3 > 1 such that for all 1 < (3 < j3 we have: 



E2 



Xo 



<c{p)\e-d'\ 



Proof. Using the expression of Zi^g^gi, and the equation (8) with the Markov 
property, it suffices to bound the quantity 



E. 



'xo,6 



(62) 



By Theorem 5 with k = 1 the function 6 — pxaiU"", V^, 0) is absolutely con- 



tinuous and we can write the quantity above as: E, 



ceo, I 
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Using first the Minkowski inequality, a change of measure and then the 
Holder inequality one finds the following bounds for this quantity: 



Je 

'-f 

Je 



El 



xo,s 



E. 



Xo,S 



pg„([/",y",g) 

p" (?7",F",s) 



ds 



p" (C/-,F",s) 



13a' 



1 

13c,' 



E. 



Xo,S 



p" (c/",y",s) 



(/3-l)/3'- 



ds 



for two conjugate exponents a' and f3' . But the first expectation in the right 
hand side above is bounded by Corollary 2 (with k = 1) for all choices of a', 
p. The second expectation can be bounded if (/3 — 1)(3' is close enough to 
zero by using (9) as in the proof of Lemma 7. This gives that (62) is smaller 
thanc\e-e'\. □ 



5.3 A technical lemma 

Lemma 10. Let (Go, . . . , Gk) be a centered Gaussian vector with invertible 
covariance matrix Gk+i and let us denote by Ck-i the covariance matrix of 
(Gi,...,Gfe_i). Then, 

X] GA^k+iXjjiGji - X Gj[Gk-i]]],Gji, (63) 

0<j,j'<k l<j,j'<k-l 
is a X^(2) random variable. 

Proof. Write the Gram-Schmidt orthonormalization procedure for the L"^ 
vectors Gi , . . . , , Go as: 



Ho 




Go 




= Pk 








_Gfc_ 



where the variables TLq., . . . , Tik are i.i.d. with standard Gaussian law and Pk 
is some triangular matrix. Then a few linear algebra shows that (63) is equal 
to Y!j=o ^1 - Y!j=i = Hl+ Hi and thus is chi-square distributed. □ 
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